\(\int (d+e x^2)^2 (a+b \sec ^{-1}(c x)) \, dx\) [82]

Optimal result

Integrand size = 18, antiderivative size = 191 \[ \int \left (d+e x^2\right )^2 \left (a+b \sec ^{-1}(c x)\right ) \, dx=-\frac {b e \left (40 c^2 d+9 e\right ) x^2 \sqrt {-1+c^2 x^2}}{120 c^3 \sqrt {c^2 x^2}}-\frac {b e^2 x^4 \sqrt {-1+c^2 x^2}}{20 c \sqrt {c^2 x^2}}+d^2 x \left (a+b \sec ^{-1}(c x)\right )+\frac {2}{3} d e x^3 \left (a+b \sec ^{-1}(c x)\right )+\frac {1}{5} e^2 x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {b \left (120 c^4 d^2+40 c^2 d e+9 e^2\right ) x \text {arctanh}\left (\frac {c x}{\sqrt {-1+c^2 x^2}}\right )}{120 c^4 \sqrt {c^2 x^2}} \]

[Out]

d^2*x*(a+b*arcsec(c*x))+2/3*d*e*x^3*(a+b*arcsec(c*x))+1/5*e^2*x^5*(a+b*arcsec(c*x))-1/120*b*(120*c^4*d^2+40*c^
2*d*e+9*e^2)*x*arctanh(c*x/(c^2*x^2-1)^(1/2))/c^4/(c^2*x^2)^(1/2)-1/120*b*e*(40*c^2*d+9*e)*x^2*(c^2*x^2-1)^(1/
2)/c^3/(c^2*x^2)^(1/2)-1/20*b*e^2*x^4*(c^2*x^2-1)^(1/2)/c/(c^2*x^2)^(1/2)

Rubi [A] (verified)

Time = 0.08 (sec) , antiderivative size = 191, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 7, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.389, Rules used = {200, 5336, 12, 1173, 396, 223, 212} \[ \int \left (d+e x^2\right )^2 \left (a+b \sec ^{-1}(c x)\right ) \, dx=d^2 x \left (a+b \sec ^{-1}(c x)\right )+\frac {2}{3} d e x^3 \left (a+b \sec ^{-1}(c x)\right )+\frac {1}{5} e^2 x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {b x \text {arctanh}\left (\frac {c x}{\sqrt {c^2 x^2-1}}\right ) \left (120 c^4 d^2+40 c^2 d e+9 e^2\right )}{120 c^4 \sqrt {c^2 x^2}}-\frac {b e^2 x^4 \sqrt {c^2 x^2-1}}{20 c \sqrt {c^2 x^2}}-\frac {b e x^2 \sqrt {c^2 x^2-1} \left (40 c^2 d+9 e\right )}{120 c^3 \sqrt {c^2 x^2}} \]

[In]

Int[(d + e*x^2)^2*(a + b*ArcSec[c*x]),x]

[Out]

-1/120*(b*e*(40*c^2*d + 9*e)*x^2*Sqrt[-1 + c^2*x^2])/(c^3*Sqrt[c^2*x^2]) - (b*e^2*x^4*Sqrt[-1 + c^2*x^2])/(20*
c*Sqrt[c^2*x^2]) + d^2*x*(a + b*ArcSec[c*x]) + (2*d*e*x^3*(a + b*ArcSec[c*x]))/3 + (e^2*x^5*(a + b*ArcSec[c*x]
))/5 - (b*(120*c^4*d^2 + 40*c^2*d*e + 9*e^2)*x*ArcTanh[(c*x)/Sqrt[-1 + c^2*x^2]])/(120*c^4*Sqrt[c^2*x^2])

Rule 12

Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] &&  !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]

Rule 200

Int[((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> Int[ExpandIntegrand[(a + b*x^n)^p, x], x] /; FreeQ[{a, b}, x]
&& IGtQ[n, 0] && IGtQ[p, 0]

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 223

Int[1/Sqrt[(a_) + (b_.)*(x_)^2], x_Symbol] :> Subst[Int[1/(1 - b*x^2), x], x, x/Sqrt[a + b*x^2]] /; FreeQ[{a,
b}, x] &&  !GtQ[a, 0]

Rule 396

Int[((a_) + (b_.)*(x_)^(n_))^(p_)*((c_) + (d_.)*(x_)^(n_)), x_Symbol] :> Simp[d*x*((a + b*x^n)^(p + 1)/(b*(n*(
p + 1) + 1))), x] - Dist[(a*d - b*c*(n*(p + 1) + 1))/(b*(n*(p + 1) + 1)), Int[(a + b*x^n)^p, x], x] /; FreeQ[{
a, b, c, d, n}, x] && NeQ[b*c - a*d, 0] && NeQ[n*(p + 1) + 1, 0]

Rule 1173

Int[((d_) + (e_.)*(x_)^2)^(q_)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_.), x_Symbol] :> Simp[c^p*x^(4*p - 1)*(
(d + e*x^2)^(q + 1)/(e*(4*p + 2*q + 1))), x] + Dist[1/(e*(4*p + 2*q + 1)), Int[(d + e*x^2)^q*ExpandToSum[e*(4*
p + 2*q + 1)*(a + b*x^2 + c*x^4)^p - d*c^p*(4*p - 1)*x^(4*p - 2) - e*c^p*(4*p + 2*q + 1)*x^(4*p), x], x], x] /
; FreeQ[{a, b, c, d, e, q}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && IGtQ[p, 0] &&  !LtQ[
q, -1]

Rule 5336

Int[((a_.) + ArcSec[(c_.)*(x_)]*(b_.))*((d_.) + (e_.)*(x_)^2)^(p_.), x_Symbol] :> With[{u = IntHide[(d + e*x^2
)^p, x]}, Dist[a + b*ArcSec[c*x], u, x] - Dist[b*c*(x/Sqrt[c^2*x^2]), Int[SimplifyIntegrand[u/(x*Sqrt[c^2*x^2
- 1]), x], x], x]] /; FreeQ[{a, b, c, d, e}, x] && (IGtQ[p, 0] || ILtQ[p + 1/2, 0])

Rubi steps \begin{align*} \text {integral}& = d^2 x \left (a+b \sec ^{-1}(c x)\right )+\frac {2}{3} d e x^3 \left (a+b \sec ^{-1}(c x)\right )+\frac {1}{5} e^2 x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {(b c x) \int \frac {15 d^2+10 d e x^2+3 e^2 x^4}{15 \sqrt {-1+c^2 x^2}} \, dx}{\sqrt {c^2 x^2}} \\ & = d^2 x \left (a+b \sec ^{-1}(c x)\right )+\frac {2}{3} d e x^3 \left (a+b \sec ^{-1}(c x)\right )+\frac {1}{5} e^2 x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {(b c x) \int \frac {15 d^2+10 d e x^2+3 e^2 x^4}{\sqrt {-1+c^2 x^2}} \, dx}{15 \sqrt {c^2 x^2}} \\ & = -\frac {b e^2 x^4 \sqrt {-1+c^2 x^2}}{20 c \sqrt {c^2 x^2}}+d^2 x \left (a+b \sec ^{-1}(c x)\right )+\frac {2}{3} d e x^3 \left (a+b \sec ^{-1}(c x)\right )+\frac {1}{5} e^2 x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {(b x) \int \frac {60 c^2 d^2+e \left (40 c^2 d+9 e\right ) x^2}{\sqrt {-1+c^2 x^2}} \, dx}{60 c \sqrt {c^2 x^2}} \\ & = -\frac {b e \left (40 c^2 d+9 e\right ) x^2 \sqrt {-1+c^2 x^2}}{120 c^3 \sqrt {c^2 x^2}}-\frac {b e^2 x^4 \sqrt {-1+c^2 x^2}}{20 c \sqrt {c^2 x^2}}+d^2 x \left (a+b \sec ^{-1}(c x)\right )+\frac {2}{3} d e x^3 \left (a+b \sec ^{-1}(c x)\right )+\frac {1}{5} e^2 x^5 \left (a+b \sec ^{-1}(c x)\right )+\frac {\left (b \left (-120 c^4 d^2-e \left (40 c^2 d+9 e\right )\right ) x\right ) \int \frac {1}{\sqrt {-1+c^2 x^2}} \, dx}{120 c^3 \sqrt {c^2 x^2}} \\ & = -\frac {b e \left (40 c^2 d+9 e\right ) x^2 \sqrt {-1+c^2 x^2}}{120 c^3 \sqrt {c^2 x^2}}-\frac {b e^2 x^4 \sqrt {-1+c^2 x^2}}{20 c \sqrt {c^2 x^2}}+d^2 x \left (a+b \sec ^{-1}(c x)\right )+\frac {2}{3} d e x^3 \left (a+b \sec ^{-1}(c x)\right )+\frac {1}{5} e^2 x^5 \left (a+b \sec ^{-1}(c x)\right )+\frac {\left (b \left (-120 c^4 d^2-e \left (40 c^2 d+9 e\right )\right ) x\right ) \text {Subst}\left (\int \frac {1}{1-c^2 x^2} \, dx,x,\frac {x}{\sqrt {-1+c^2 x^2}}\right )}{120 c^3 \sqrt {c^2 x^2}} \\ & = -\frac {b e \left (40 c^2 d+9 e\right ) x^2 \sqrt {-1+c^2 x^2}}{120 c^3 \sqrt {c^2 x^2}}-\frac {b e^2 x^4 \sqrt {-1+c^2 x^2}}{20 c \sqrt {c^2 x^2}}+d^2 x \left (a+b \sec ^{-1}(c x)\right )+\frac {2}{3} d e x^3 \left (a+b \sec ^{-1}(c x)\right )+\frac {1}{5} e^2 x^5 \left (a+b \sec ^{-1}(c x)\right )-\frac {b \left (120 c^4 d^2+40 c^2 d e+9 e^2\right ) x \text {arctanh}\left (\frac {c x}{\sqrt {-1+c^2 x^2}}\right )}{120 c^4 \sqrt {c^2 x^2}} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.14 (sec) , antiderivative size = 153, normalized size of antiderivative = 0.80 \[ \int \left (d+e x^2\right )^2 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\frac {c^2 x \left (8 a c^3 \left (15 d^2+10 d e x^2+3 e^2 x^4\right )-b e \sqrt {1-\frac {1}{c^2 x^2}} x \left (9 e+c^2 \left (40 d+6 e x^2\right )\right )\right )+8 b c^5 x \left (15 d^2+10 d e x^2+3 e^2 x^4\right ) \sec ^{-1}(c x)-b \left (120 c^4 d^2+40 c^2 d e+9 e^2\right ) \log \left (\left (1+\sqrt {1-\frac {1}{c^2 x^2}}\right ) x\right )}{120 c^5} \]

[In]

Integrate[(d + e*x^2)^2*(a + b*ArcSec[c*x]),x]

[Out]

(c^2*x*(8*a*c^3*(15*d^2 + 10*d*e*x^2 + 3*e^2*x^4) - b*e*Sqrt[1 - 1/(c^2*x^2)]*x*(9*e + c^2*(40*d + 6*e*x^2)))
+ 8*b*c^5*x*(15*d^2 + 10*d*e*x^2 + 3*e^2*x^4)*ArcSec[c*x] - b*(120*c^4*d^2 + 40*c^2*d*e + 9*e^2)*Log[(1 + Sqrt
[1 - 1/(c^2*x^2)])*x])/(120*c^5)

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(339\) vs. \(2(169)=338\).

Time = 0.42 (sec) , antiderivative size = 340, normalized size of antiderivative = 1.78

method	result	size
parts	\(a \left (\frac {1}{5} e^{2} x^{5}+\frac {2}{3} d e \,x^{3}+d^{2} x \right )+\frac {b \,\operatorname {arcsec}\left (c x \right ) e^{2} x^{5}}{5}+\frac {2 b \,\operatorname {arcsec}\left (c x \right ) d e \,x^{3}}{3}+b \,\operatorname {arcsec}\left (c x \right ) d^{2} x -\frac {b \left (c^{2} x^{2}-1\right ) x^{2} e^{2}}{20 c^{3} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {b \left (c^{2} x^{2}-1\right ) d e}{3 c^{3} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {b \sqrt {c^{2} x^{2}-1}\, d^{2} \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{c^{2} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, x}-\frac {3 b \left (c^{2} x^{2}-1\right ) e^{2}}{40 c^{5} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {b \sqrt {c^{2} x^{2}-1}\, d e \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{3 c^{4} x \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {3 b \sqrt {c^{2} x^{2}-1}\, e^{2} \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{40 c^{6} x \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}\)	\(340\)
derivativedivides	\(\frac {\frac {a \left (d^{2} c^{5} x +\frac {2}{3} d \,c^{5} e \,x^{3}+\frac {1}{5} e^{2} c^{5} x^{5}\right )}{c^{4}}+b \,\operatorname {arcsec}\left (c x \right ) d^{2} c x +\frac {2 b c \,\operatorname {arcsec}\left (c x \right ) d e \,x^{3}}{3}+\frac {b c \,\operatorname {arcsec}\left (c x \right ) e^{2} x^{5}}{5}-\frac {b \sqrt {c^{2} x^{2}-1}\, d^{2} \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{\sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, c x}-\frac {b \left (c^{2} x^{2}-1\right ) d e}{3 c^{2} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {b \left (c^{2} x^{2}-1\right ) x^{2} e^{2}}{20 c^{2} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {b \sqrt {c^{2} x^{2}-1}\, d e \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{3 c^{3} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, x}-\frac {3 b \left (c^{2} x^{2}-1\right ) e^{2}}{40 c^{4} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {3 b \sqrt {c^{2} x^{2}-1}\, e^{2} \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{40 c^{5} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, x}}{c}\)	\(359\)
default	\(\frac {\frac {a \left (d^{2} c^{5} x +\frac {2}{3} d \,c^{5} e \,x^{3}+\frac {1}{5} e^{2} c^{5} x^{5}\right )}{c^{4}}+b \,\operatorname {arcsec}\left (c x \right ) d^{2} c x +\frac {2 b c \,\operatorname {arcsec}\left (c x \right ) d e \,x^{3}}{3}+\frac {b c \,\operatorname {arcsec}\left (c x \right ) e^{2} x^{5}}{5}-\frac {b \sqrt {c^{2} x^{2}-1}\, d^{2} \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{\sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, c x}-\frac {b \left (c^{2} x^{2}-1\right ) d e}{3 c^{2} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {b \left (c^{2} x^{2}-1\right ) x^{2} e^{2}}{20 c^{2} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {b \sqrt {c^{2} x^{2}-1}\, d e \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{3 c^{3} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, x}-\frac {3 b \left (c^{2} x^{2}-1\right ) e^{2}}{40 c^{4} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}}-\frac {3 b \sqrt {c^{2} x^{2}-1}\, e^{2} \ln \left (c x +\sqrt {c^{2} x^{2}-1}\right )}{40 c^{5} \sqrt {\frac {c^{2} x^{2}-1}{c^{2} x^{2}}}\, x}}{c}\)	\(359\)

[In]

int((e*x^2+d)^2*(a+b*arcsec(c*x)),x,method=_RETURNVERBOSE)

[Out]

a*(1/5*e^2*x^5+2/3*d*e*x^3+d^2*x)+1/5*b*arcsec(c*x)*e^2*x^5+2/3*b*arcsec(c*x)*d*e*x^3+b*arcsec(c*x)*d^2*x-1/20
*b/c^3*(c^2*x^2-1)*x^2/((c^2*x^2-1)/c^2/x^2)^(1/2)*e^2-1/3*b/c^3*(c^2*x^2-1)/((c^2*x^2-1)/c^2/x^2)^(1/2)*d*e-b
/c^2*(c^2*x^2-1)^(1/2)/((c^2*x^2-1)/c^2/x^2)^(1/2)/x*d^2*ln(c*x+(c^2*x^2-1)^(1/2))-3/40*b/c^5*(c^2*x^2-1)/((c^
2*x^2-1)/c^2/x^2)^(1/2)*e^2-1/3*b/c^4*(c^2*x^2-1)^(1/2)/x/((c^2*x^2-1)/c^2/x^2)^(1/2)*d*e*ln(c*x+(c^2*x^2-1)^(
1/2))-3/40*b/c^6*(c^2*x^2-1)^(1/2)/x/((c^2*x^2-1)/c^2/x^2)^(1/2)*e^2*ln(c*x+(c^2*x^2-1)^(1/2))

Fricas [A] (verification not implemented)

none

Time = 0.39 (sec) , antiderivative size = 237, normalized size of antiderivative = 1.24 \[ \int \left (d+e x^2\right )^2 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\frac {24 \, a c^{5} e^{2} x^{5} + 80 \, a c^{5} d e x^{3} + 120 \, a c^{5} d^{2} x + 8 \, {\left (3 \, b c^{5} e^{2} x^{5} + 10 \, b c^{5} d e x^{3} + 15 \, b c^{5} d^{2} x - 15 \, b c^{5} d^{2} - 10 \, b c^{5} d e - 3 \, b c^{5} e^{2}\right )} \operatorname {arcsec}\left (c x\right ) + 16 \, {\left (15 \, b c^{5} d^{2} + 10 \, b c^{5} d e + 3 \, b c^{5} e^{2}\right )} \arctan \left (-c x + \sqrt {c^{2} x^{2} - 1}\right ) + {\left (120 \, b c^{4} d^{2} + 40 \, b c^{2} d e + 9 \, b e^{2}\right )} \log \left (-c x + \sqrt {c^{2} x^{2} - 1}\right ) - {\left (6 \, b c^{3} e^{2} x^{3} + {\left (40 \, b c^{3} d e + 9 \, b c e^{2}\right )} x\right )} \sqrt {c^{2} x^{2} - 1}}{120 \, c^{5}} \]

[In]

integrate((e*x^2+d)^2*(a+b*arcsec(c*x)),x, algorithm="fricas")

[Out]

1/120*(24*a*c^5*e^2*x^5 + 80*a*c^5*d*e*x^3 + 120*a*c^5*d^2*x + 8*(3*b*c^5*e^2*x^5 + 10*b*c^5*d*e*x^3 + 15*b*c^
5*d^2*x - 15*b*c^5*d^2 - 10*b*c^5*d*e - 3*b*c^5*e^2)*arcsec(c*x) + 16*(15*b*c^5*d^2 + 10*b*c^5*d*e + 3*b*c^5*e
^2)*arctan(-c*x + sqrt(c^2*x^2 - 1)) + (120*b*c^4*d^2 + 40*b*c^2*d*e + 9*b*e^2)*log(-c*x + sqrt(c^2*x^2 - 1))
- (6*b*c^3*e^2*x^3 + (40*b*c^3*d*e + 9*b*c*e^2)*x)*sqrt(c^2*x^2 - 1))/c^5

Sympy [A] (verification not implemented)

Time = 5.83 (sec) , antiderivative size = 355, normalized size of antiderivative = 1.86 \[ \int \left (d+e x^2\right )^2 \left (a+b \sec ^{-1}(c x)\right ) \, dx=a d^{2} x + \frac {2 a d e x^{3}}{3} + \frac {a e^{2} x^{5}}{5} + b d^{2} x \operatorname {asec}{\left (c x \right )} + \frac {2 b d e x^{3} \operatorname {asec}{\left (c x \right )}}{3} + \frac {b e^{2} x^{5} \operatorname {asec}{\left (c x \right )}}{5} - \frac {b d^{2} \left (\begin {cases} \operatorname {acosh}{\left (c x \right )} & \text {for}\: \left |{c^{2} x^{2}}\right | > 1 \\- i \operatorname {asin}{\left (c x \right )} & \text {otherwise} \end {cases}\right )}{c} - \frac {2 b d e \left (\begin {cases} \frac {x \sqrt {c^{2} x^{2} - 1}}{2 c} + \frac {\operatorname {acosh}{\left (c x \right )}}{2 c^{2}} & \text {for}\: \left |{c^{2} x^{2}}\right | > 1 \\- \frac {i c x^{3}}{2 \sqrt {- c^{2} x^{2} + 1}} + \frac {i x}{2 c \sqrt {- c^{2} x^{2} + 1}} - \frac {i \operatorname {asin}{\left (c x \right )}}{2 c^{2}} & \text {otherwise} \end {cases}\right )}{3 c} - \frac {b e^{2} \left (\begin {cases} \frac {c x^{5}}{4 \sqrt {c^{2} x^{2} - 1}} + \frac {x^{3}}{8 c \sqrt {c^{2} x^{2} - 1}} - \frac {3 x}{8 c^{3} \sqrt {c^{2} x^{2} - 1}} + \frac {3 \operatorname {acosh}{\left (c x \right )}}{8 c^{4}} & \text {for}\: \left |{c^{2} x^{2}}\right | > 1 \\- \frac {i c x^{5}}{4 \sqrt {- c^{2} x^{2} + 1}} - \frac {i x^{3}}{8 c \sqrt {- c^{2} x^{2} + 1}} + \frac {3 i x}{8 c^{3} \sqrt {- c^{2} x^{2} + 1}} - \frac {3 i \operatorname {asin}{\left (c x \right )}}{8 c^{4}} & \text {otherwise} \end {cases}\right )}{5 c} \]

[In]

integrate((e*x**2+d)**2*(a+b*asec(c*x)),x)

[Out]

a*d**2*x + 2*a*d*e*x**3/3 + a*e**2*x**5/5 + b*d**2*x*asec(c*x) + 2*b*d*e*x**3*asec(c*x)/3 + b*e**2*x**5*asec(c
*x)/5 - b*d**2*Piecewise((acosh(c*x), Abs(c**2*x**2) > 1), (-I*asin(c*x), True))/c - 2*b*d*e*Piecewise((x*sqrt
(c**2*x**2 - 1)/(2*c) + acosh(c*x)/(2*c**2), Abs(c**2*x**2) > 1), (-I*c*x**3/(2*sqrt(-c**2*x**2 + 1)) + I*x/(2
*c*sqrt(-c**2*x**2 + 1)) - I*asin(c*x)/(2*c**2), True))/(3*c) - b*e**2*Piecewise((c*x**5/(4*sqrt(c**2*x**2 - 1
)) + x**3/(8*c*sqrt(c**2*x**2 - 1)) - 3*x/(8*c**3*sqrt(c**2*x**2 - 1)) + 3*acosh(c*x)/(8*c**4), Abs(c**2*x**2)
> 1), (-I*c*x**5/(4*sqrt(-c**2*x**2 + 1)) - I*x**3/(8*c*sqrt(-c**2*x**2 + 1)) + 3*I*x/(8*c**3*sqrt(-c**2*x**2
+ 1)) - 3*I*asin(c*x)/(8*c**4), True))/(5*c)

Maxima [A] (verification not implemented)

none

Time = 0.23 (sec) , antiderivative size = 296, normalized size of antiderivative = 1.55 \[ \int \left (d+e x^2\right )^2 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\frac {1}{5} \, a e^{2} x^{5} + \frac {2}{3} \, a d e x^{3} + \frac {1}{6} \, {\left (4 \, x^{3} \operatorname {arcsec}\left (c x\right ) - \frac {\frac {2 \, \sqrt {-\frac {1}{c^{2} x^{2}} + 1}}{c^{2} {\left (\frac {1}{c^{2} x^{2}} - 1\right )} + c^{2}} + \frac {\log \left (\sqrt {-\frac {1}{c^{2} x^{2}} + 1} + 1\right )}{c^{2}} - \frac {\log \left (\sqrt {-\frac {1}{c^{2} x^{2}} + 1} - 1\right )}{c^{2}}}{c}\right )} b d e + \frac {1}{80} \, {\left (16 \, x^{5} \operatorname {arcsec}\left (c x\right ) + \frac {\frac {2 \, {\left (3 \, {\left (-\frac {1}{c^{2} x^{2}} + 1\right )}^{\frac {3}{2}} - 5 \, \sqrt {-\frac {1}{c^{2} x^{2}} + 1}\right )}}{c^{4} {\left (\frac {1}{c^{2} x^{2}} - 1\right )}^{2} + 2 \, c^{4} {\left (\frac {1}{c^{2} x^{2}} - 1\right )} + c^{4}} - \frac {3 \, \log \left (\sqrt {-\frac {1}{c^{2} x^{2}} + 1} + 1\right )}{c^{4}} + \frac {3 \, \log \left (\sqrt {-\frac {1}{c^{2} x^{2}} + 1} - 1\right )}{c^{4}}}{c}\right )} b e^{2} + a d^{2} x + \frac {{\left (2 \, c x \operatorname {arcsec}\left (c x\right ) - \log \left (\sqrt {-\frac {1}{c^{2} x^{2}} + 1} + 1\right ) + \log \left (-\sqrt {-\frac {1}{c^{2} x^{2}} + 1} + 1\right )\right )} b d^{2}}{2 \, c} \]

[In]

integrate((e*x^2+d)^2*(a+b*arcsec(c*x)),x, algorithm="maxima")

[Out]

1/5*a*e^2*x^5 + 2/3*a*d*e*x^3 + 1/6*(4*x^3*arcsec(c*x) - (2*sqrt(-1/(c^2*x^2) + 1)/(c^2*(1/(c^2*x^2) - 1) + c^
2) + log(sqrt(-1/(c^2*x^2) + 1) + 1)/c^2 - log(sqrt(-1/(c^2*x^2) + 1) - 1)/c^2)/c)*b*d*e + 1/80*(16*x^5*arcsec
(c*x) + (2*(3*(-1/(c^2*x^2) + 1)^(3/2) - 5*sqrt(-1/(c^2*x^2) + 1))/(c^4*(1/(c^2*x^2) - 1)^2 + 2*c^4*(1/(c^2*x^
2) - 1) + c^4) - 3*log(sqrt(-1/(c^2*x^2) + 1) + 1)/c^4 + 3*log(sqrt(-1/(c^2*x^2) + 1) - 1)/c^4)/c)*b*e^2 + a*d
^2*x + 1/2*(2*c*x*arcsec(c*x) - log(sqrt(-1/(c^2*x^2) + 1) + 1) + log(-sqrt(-1/(c^2*x^2) + 1) + 1))*b*d^2/c

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 14166 vs. \(2 (169) = 338\).

Time = 4.60 (sec) , antiderivative size = 14166, normalized size of antiderivative = 74.17 \[ \int \left (d+e x^2\right )^2 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\text {Too large to display} \]

[In]

integrate((e*x^2+d)^2*(a+b*arcsec(c*x)),x, algorithm="giac")

[Out]

1/120*(120*b*c^4*d^2*arccos(1/(c*x))/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)
^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 +
c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10) - 120*b*c^4*d^2*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/(c^6
+ 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) -
1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)
+ 120*b*c^4*d^2*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2
+ 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2)
- 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10) + 120*a*c^4*d^2/(c^6 + 5*c^6*(1/(c^2*x^2)
- 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6
+ 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10) + 360*b*c^4*d^2*(1/(c
^2*x^2) - 1)*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(
c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/
(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) - 600*b*c^4*d^2*(1/(c^2*x^2) - 1)*log(abs(sqrt(-1/(c^2*x^2
) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) +
1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x
^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) + 600*b*c^4*d^2*(1/(c^2*x^2) - 1)*log(abs(sqrt(-1/(c^2*x^2) + 1)
- 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4
+ 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) -
1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) + 360*a*c^4*d^2*(1/(c^2*x^2) - 1)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(
c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*
(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) + 80*b*c^2*d*
e*arccos(1/(c*x))/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4
+ 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) -
1)^5/(1/(c*x) + 1)^10) + 240*b*c^4*d^2*(1/(c^2*x^2) - 1)^2*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/
(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6
*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) - 40*b*c^2*d
*e*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(
c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c
*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10) - 1200*b*c^4*d^2*(1/(c^2*x^2) - 1)^2*log(abs(sqrt(-1/(c
^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(
c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/
(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) + 40*b*c^2*d*e*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) -
1))/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c
^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x)
+ 1)^10) + 1200*b*c^4*d^2*(1/(c^2*x^2) - 1)^2*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/
(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c
*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) +
1)^4) + 80*a*c^2*d*e/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)
^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2)
- 1)^5/(1/(c*x) + 1)^10) + 240*a*c^4*d^2*(1/(c^2*x^2) - 1)^2/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2
+ 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2)
- 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) - 80*b*c^2*d*e*(1/(c^2*x^
2) - 1)*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x)
+ 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*
x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) - 240*b*c^4*d^2*(1/(c^2*x^2) - 1)^3*arccos(1/(c*x))/((c^6 + 5*c
^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3
/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c
*x) + 1)^6) - 200*b*c^2*d*e*(1/(c^2*x^2) - 1)*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/
(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c
*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) +
1)^2) - 1200*b*c^4*d^2*(1/(c^2*x^2) - 1)^3*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^
2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x)
+ 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^
6) + 200*b*c^2*d*e*(1/(c^2*x^2) - 1)*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2)
- 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^
6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) + 1
200*b*c^4*d^2*(1/(c^2*x^2) - 1)^3*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) -
1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 +
5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) - 80*b
*c^2*d*e*sqrt(-1/(c^2*x^2) + 1)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(
1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*
(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)) - 80*a*c^2*d*e*(1/(c^2*x^2) - 1)/((c^6 + 5*c^6*(1/(c^2*x^
2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1
)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) -
240*a*c^4*d^2*(1/(c^2*x^2) - 1)^3/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^
2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c
^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) + 24*b*e^2*arccos(1/(c*x))/(c^6 + 5*c^6*(1/(c^2*x^2)
- 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^
6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10) - 160*b*c^2*d*e*(1/(
c^2*x^2) - 1)^2*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(
1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*
(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) - 360*b*c^4*d^2*(1/(c^2*x^2) - 1)^4*arccos(1/(c*x))/((c
^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2
) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^1
0)*(1/(c*x) + 1)^8) - 9*b*e^2*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1
/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^
6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10) - 400*b*c^2*d*e*(1/(c^2*x^2)
- 1)^2*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^
6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4
/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) - 600*b*c^4*d^2*(1/(c^2*x^2) - 1
)^4*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1
/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/
(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^8) + 9*b*e^2*log(abs(sqrt(-1/(c^2*x^2)
+ 1) - 1/(c*x) - 1))/(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)
^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2)
- 1)^5/(1/(c*x) + 1)^10) + 400*b*c^2*d*e*(1/(c^2*x^2) - 1)^2*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/(
(c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x
^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)
^10)*(1/(c*x) + 1)^4) + 600*b*c^4*d^2*(1/(c^2*x^2) - 1)^4*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6
+ 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2)
- 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)
*(1/(c*x) + 1)^8) + 160*b*c^2*d*e*(-1/(c^2*x^2) + 1)^(3/2)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 1
0*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) -
1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^3) + 24*a*e^2/(c^6 + 5*c^6*(1/(
c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*
x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10) - 160*a*c^2*
d*e*(1/(c^2*x^2) - 1)^2/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x)
+ 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*
x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) - 360*a*c^4*d^2*(1/(c^2*x^2) - 1)^4/((c^6 + 5*c^6*(1/(c^2*x^2)
- 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6
+ 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^8) - 12
0*b*e^2*(1/(c^2*x^2) - 1)*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2
) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1
)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) + 160*b*c^2*d*e*(1/(c^2*x^2) - 1)^3*arccos(1/
(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(
1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c
*x) + 1)^10)*(1/(c*x) + 1)^6) - 120*b*c^4*d^2*(1/(c^2*x^2) - 1)^5*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) -
1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6
+ 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^10) - 45
*b*e^2*(1/(c^2*x^2) - 1)*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*
x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1
/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) - 400*b*c^2*d*e
*(1/(c^2*x^2) - 1)^3*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) +
1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^
2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) - 120*b*c^4*d^2*(1/
(c^2*x^2) - 1)^5*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^
2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^
2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^10) + 45*b*e^2*(1/(c^2*x^2
) - 1)*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6
*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/
(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) + 400*b*c^2*d*e*(1/(c^2*x^2) - 1)
^3*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/
(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(
c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) + 120*b*c^4*d^2*(1/(c^2*x^2) - 1)^5*l
og(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2
*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x)
+ 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^10) - 30*b*e^2*sqrt(-1/(c^2*x^2) + 1)/((c^6
+ 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) -
1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*
(1/(c*x) + 1)) - 120*a*e^2*(1/(c^2*x^2) - 1)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*
x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x)
+ 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^2) + 160*a*c^2*d*e*(1/(c^2*x^2) - 1)^3/((c^6
+ 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) -
1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*
(1/(c*x) + 1)^6) - 120*a*c^4*d^2*(1/(c^2*x^2) - 1)^5/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*
(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(
1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^10) + 240*b*e^2*(1/(c^2*x^2) - 1)^2*a
rccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 +
10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)
^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) + 80*b*c^2*d*e*(1/(c^2*x^2) - 1)^4*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2
*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x)
+ 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^8
) - 90*b*e^2*(1/(c^2*x^2) - 1)^2*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1
)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 +
5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) - 200*b
*c^2*d*e*(1/(c^2*x^2) - 1)^4*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1
/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^
6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^8) + 90*b*e^2*
(1/(c^2*x^2) - 1)^2*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) +
1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2
*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) + 200*b*c^2*d*e*(1/(
c^2*x^2) - 1)^4*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2
+ 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2
) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^8) + 160*b*c^2*d*e*(1/(c^2*
x^2) - 1)^3*sqrt(-1/(c^2*x^2) + 1)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^
2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c
^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^7) + 12*b*e^2*(-1/(c^2*x^2) + 1)^(3/2)/((c^6 + 5*c^6*(1
/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(
c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) +
1)^3) + 240*a*e^2*(1/(c^2*x^2) - 1)^2/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) -
1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8
+ c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^4) + 80*a*c^2*d*e*(1/(c^2*x^2) - 1)^4/((c^6 + 5*c^6
*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(
1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x
) + 1)^8) - 240*b*e^2*(1/(c^2*x^2) - 1)^3*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10
*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1
)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) - 80*b*c^2*d*e*(1/(c^2*x^2) -
1)^5*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) +
1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^
2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^10) - 90*b*e^2*(1/(c^2*x^2) - 1)^3*log(abs(sqrt(-1/(c^2*x^2) + 1) +
1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 1
0*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^
5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) - 40*b*c^2*d*e*(1/(c^2*x^2) - 1)^5*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*
x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6
*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/
(c*x) + 1)^10)*(1/(c*x) + 1)^10) + 90*b*e^2*(1/(c^2*x^2) - 1)^3*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))
/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2
*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) +
1)^10)*(1/(c*x) + 1)^6) + 40*b*c^2*d*e*(1/(c^2*x^2) - 1)^5*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^
6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2)
- 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10
)*(1/(c*x) + 1)^10) + 80*b*c^2*d*e*(1/(c^2*x^2) - 1)^4*sqrt(-1/(c^2*x^2) + 1)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/
(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*
c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^9) - 240*a*e
^2*(1/(c^2*x^2) - 1)^3/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) +
1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x
^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^6) - 80*a*c^2*d*e*(1/(c^2*x^2) - 1)^5/((c^6 + 5*c^6*(1/(c^2*x^2) -
1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 +
5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^10) + 120
*b*e^2*(1/(c^2*x^2) - 1)^4*arccos(1/(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^
2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) +
1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^8) - 45*b*e^2*(1/(c^2*x^2) - 1)^4*log(abs(sqrt(
-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2
/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^
6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^8) + 45*b*e^2*(1/(c^2*x^2) - 1)^4*log(abs(sqrt(-1/(c^2*x
^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x)
+ 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2
*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^8) + 12*b*e^2*(1/(c^2*x^2) - 1)^3*sqrt(-1/(c^2*x^2) + 1)/((c^6 +
5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1
)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1
/(c*x) + 1)^7) + 120*a*e^2*(1/(c^2*x^2) - 1)^4/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^
2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x
) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^8) - 24*b*e^2*(1/(c^2*x^2) - 1)^5*arccos(1/
(c*x))/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(
1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c
*x) + 1)^10)*(1/(c*x) + 1)^10) - 9*b*e^2*(1/(c^2*x^2) - 1)^5*log(abs(sqrt(-1/(c^2*x^2) + 1) + 1/(c*x) + 1))/((
c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^
2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^
10)*(1/(c*x) + 1)^10) + 9*b*e^2*(1/(c^2*x^2) - 1)^5*log(abs(sqrt(-1/(c^2*x^2) + 1) - 1/(c*x) - 1))/((c^6 + 5*c
^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3
/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c
*x) + 1)^10) + 30*b*e^2*(1/(c^2*x^2) - 1)^4*sqrt(-1/(c^2*x^2) + 1)/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) +
1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2
*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(1/(c*x) + 1)^10)*(1/(c*x) + 1)^9) - 24*a*e^2*(1/(c^2*x
^2) - 1)^5/((c^6 + 5*c^6*(1/(c^2*x^2) - 1)/(1/(c*x) + 1)^2 + 10*c^6*(1/(c^2*x^2) - 1)^2/(1/(c*x) + 1)^4 + 10*c
^6*(1/(c^2*x^2) - 1)^3/(1/(c*x) + 1)^6 + 5*c^6*(1/(c^2*x^2) - 1)^4/(1/(c*x) + 1)^8 + c^6*(1/(c^2*x^2) - 1)^5/(
1/(c*x) + 1)^10)*(1/(c*x) + 1)^10))*c

Mupad [F(-1)]

Timed out. \[ \int \left (d+e x^2\right )^2 \left (a+b \sec ^{-1}(c x)\right ) \, dx=\int {\left (e\,x^2+d\right )}^2\,\left (a+b\,\mathrm {acos}\left (\frac {1}{c\,x}\right )\right ) \,d x \]

[In]

int((d + e*x^2)^2*(a + b*acos(1/(c*x))),x)

[Out]

int((d + e*x^2)^2*(a + b*acos(1/(c*x))), x)